CPA for charge ordering in the extended Hubbard model 
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We study charge ordering in the extended Hubbard model 
with both on-site and nearest neighbor Coulomb repulsion 
(U and V, respectively) within the Coherent potential ap- 
proximation (CPA). The phase boundary between the homo- 
geneous and charge ordered phase for the square lattice is 
obtained for different values of U. It is shown that at quar- 
ter filling for all values of U the charge ordering exists only 
if the inter-site Coulomb repulsion V exceeds certain critical 
value which is of the order of the kinetic energy t. At finite 
temperature a reentrant transition is found in some region of 
V. 



PACS 71.10.Fd, 71.27.+a 



I. INTRODUCTION 

The problem of charge ordering (CO) attracted the 
attention of physicists already at the end of the thir- 
ties. In the low density limit, as was first proposed 
by Wignera, the electrons crystallize in form of a lat- 
tice in order to keep the Coulomb repulsion as small 
as possible. Such a Wigner lattice is experimentally 
realized in a GaAs/AlGaAs heterostructurea. The CO 
may also occur at higher electron concentration if the in- 
teraction of electrons with spin degrees of freedom or 
with phonons drastically reduces the kinetic energyu. 
Recently the CO is extensively observed in real ma- 
terials at high densities: hole ordering in rare-earth 
pnictides like Yb4As3Q, charge ordering in the uncon- 
ventional spin-Peierls material a' — NaV^OsU and in 
colossal magnetoresistance compounds, for examnk in 
R2-2 X Ai +2x Mn 2 07(R = La, Pr; A = Ca, Sr;x > 0.5)0. 

One of the simplest models of interacting electrons 
that allows for charge ordering is the extended Hubbard 
model (EHM). This model has been intensively studied 
both in low dimensions and in the limit of infinite di- 
mension, usually at half or at quarter filling. A vari- 
ety of techniques, such as Hartree-Fock approximationQ, 
perturbation theorycl, the dynamical mean field theory 
(DMFT)la, slave boson approachEJ, as well as numerical 
methods like Quaittum Monte Carlo simulationO and 
Lanczos technique^ have been employed. Obviously, 
each of these approaches is able to describe in a proper 
way only one of the relevant limits of the model and 



despite many publications devoted to the CO in solids 
the physical picture of the phenomenon is far from being 
clear. 

Recently a melting of the charge ordered state 
on decreasing the temperature has been found in 
the doped compounds Pro.6s(Cao.7Sro.;Av9fiMn03E3 and 
La 2 -2xSri + 2xMn 2 7 (0.47 < x < 0.62)BE1 A reentrant 
transition at quarter filling has been obtained theoretic 
cally using EHM both with electron-phonon interaction!!^ 
and without clcctron-phonon interactiono ' . The present 
paper is devoted to study of the boundary between the 
charge ordered and disordered phase for different regimes 
of the temperature T, the Coulomb interactions U, V and 
the band filling n. A simple but physically meaning- 
ful approximation allowing to solve this problem is the 
CPA. This self-consistent approximation is recognized as 
the best single-site approximation for the spectral prop- 
erties of disordered systems. Originally, the alloy-analog 
approximation has been formulated as an approximation 
scheme for the Hubbard modellLZl. To solve the alloy prob- 
lem the CPA is used as a second step. The CPA was 
also applied to intermediate valence and heavy fermion 
systemsEEI. In the present work this approximation is 
used for the first time to treat the charge ordering in the 
EHM. 



II. MODEL AND FORMALISM 



We consider the following Hamiltonian for the EHM: 



<ij>cr 

U^n^riii + V ^2 n i n j 



(1) 



where Cj CT (c^.) annihilates (creates) an electron with spin 
a at site i, rii a = cf^Cia and rii = n^+nii. < ij > denotes 
nearest neighbors, t is the hopping parameter, U and V 
are on-site and inter-site Coulomb repulsion, respectively. 
We divide the hypercubic lattice in two sublattices such 
that points on one sublattice have only points of the other 
sublattice as nearest neighbors. The sublattice is denoted 
by subindex A or B: = aj cr (6j cr ) if i G A (i 6 B). 
Performing a mean-field decoupling of the V term, we 
get 
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H = zVn B a+ a a ia + U}n^nji + 

Y zVn A b+ a b ]a + U ^njinji (2) 

j&B,a j&B 
<ij><y 

where z is the number of nearest neighbors, n A / B is 
the averaged electron occupation number in the A/B- 
sublattice, N is the number of sites in the lattice. 
In the alloy-analog approach the many-body Hamilto- 
nian (|J) is replaced by a one-particle Hamiltonian with 
disorder which is of the form 

H = Y E Aa a\ a a ia + Y E Br7 b+ a b jtT 



i£A,cr 



+t Y ( a t b l° + b t a ^ - \zNVn A n B (3) 



<ij>a 



where 



E 



(zVn B / A with probability 1 — n A /B,- 



A l B - a yzVn B / A + U with probability n A / B _ a 

(4) 

In the following we assume spin-independent expectation 
values in (^), i.e. we consider only nonmagnetic solution: 
n a f = n a i = \n a \ (a = A, B). The Green function G 
corresponding to the Hamiltonian (^) has to be averaged 
over all possible configurations of the random potential 
which can be considered to be due to alloy constituents. 
The averaging cannot be performed exactly. To solve 
the alloy problem the CPA is used. The averaged Green 
function G is obtained from an effective Hamiltonian con- 
taining a self-energy T, A / b (uj) for the A/B-sublattice: 

H e ff = £UM Y a t<j a ™ + H b % b J° 

+t Y ( a i<r b i° + b t* ai °) ( 5 ) 
<ij>0 

In momentum space the averaged Green functions 
G A / B (k,uj) for the A/B-sublattice are of the form 



G A {k,Lo) = [oj - Y, a (uj) 



G B (k,uj) = I u> - T, B (u) 



ft 



UJ - 

uj - T, a (uj) 



(6) 
(7) 



where tj: is the Fourier transform of the hopping matrix 
element (the wavevector k is for sublattice A (or B)), 
^a/b( w ) are to be determined latter. 
The density of states (DOS) for free electrons with 
the band dispersion tz is replaced by the semi-elliptical 



Po{u) — -^n^W 2 — uj 2 (we set W = 1 as the unit for the 
energy scale). Note that this model DOS is often used 
as an additional approximation in combination with the 
CPA. It was noted in Refo that for the Bethe lattice 
with 3 < z < 6 this approximation is a good one, at 
least in a qualitative sense. The averaged Green func- 
tions G A / B (uj) for the A/B-sublattice then take the form 

1 - 2 / 

GU(w) = — ^GA(fc,w) = ^[uj-E b (uj) 



(uj - E B (w)) 



W 2 

2 LU-Y> B (uj) 



1/2 



(8) 



And G B (ui) is obtained by replacing A «-> B. A scat- 
tering matrix T is introduced for each configuration via 



G = G + GTG 



(9) 



The CPA demands that the scattering matrix vanishes 
on average: T = 0. This yields expression for T, A i b (uj) 
of the form 



Ea(w) = E A - (zVn B - Y, A (u)) ■ 

G A (Lu)(zVn B + U -Z a {lu)) 



(10) 



where E A = zVn B + \\Jn A . Again T, b (uj) and E B are 
obtained by replacing A <-> B. 

For arbitrary size of electron density n we make the fol- 
lowing ansatz: 

n A / B =n±x; G a/b (uj) = G(±x, uj) 



Eliminating Ea(w),Eb(u) from QSJ) and ( |10[ ) leads to an 
equation for G(±x,ui): 

G(-x,u) 1 .,, v U, 

uj — = zV(n - x) + — (n + x) - 

4 G(x,ui) 2 



zV(n — x) — uj + 



G(-x,uj) 



1 



zVln - x) + U - uj + 



4 G(x,uj) 
G(~x,uj) 



1 



G(x,u) 



(11) 

G{x,uj) 



Setting x — in eq.(|ll|) and shifting the one-electron 
energy by zVn + ^ we reproduce the CPA equationJbr 
the Green function obtained by Velicky et al in Ref.E!. 
If fi is the chemical potential of electrons, then at the 
temperature T one has 



k,n 

2 r+oc 

IT 



<jw/(a-')3 ! G(n Q , — n, uj) 



(12) 



where uj n = (2n + l)7rT is the Matsubara frequencies and 
f(uj) = (1 + exp(w — ^yr) -1 is the Fermi function. 
The pair of equations (J12[) must now be solved with the 
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FIG. 1. (V-n)-phase diagram for the 2D extended Hub- 
bard model (W = 1, T = 0) for different values of 
U : (7 = 0, 0.5, 1.5 and U = oo corresponding to the dotted, 
long-dashed, dashed and solid lines, respectively. 



constraint tia + ns = 2ri for tia,tib and \i. For small 
enough V the solution of ([l2|) is the charge disordered 
state with tia = tib- But if V is sufficiently large it 
may also be possible to find the CO solution for which 
tia 7^ nB- One finds that the condition for the onset of 
charge ordering is equivalent to tia = ns = n being a 
double solution of (|l2|). This condition is expressed as 



where G"(0,w) 
of (|ll|) when x 



2 


r+oo 




y d 


n 


J — oo 


2 


r+oc 




y d 


7T 


J — oo 








dx 



duf(uj)^G(0,Lu) 



(13) 
(14) 



c =o an d G(0,uj) is a solution 



— 0. The latter is a cubic equation 
for G(0, oj) and the correct root must be identified from 
the physical condition to yield a non-negative density of 
states. It is also easy to obtain G'(0,ui) from eq. 
So, for fixed temperature T, on-site Coulomb repultion 
U and band filling n we have the closed system of equa- 



tions (13)-(Q_4J) for the critical value V and the chemical 



potential /i within the framework of the CPA. 



III. NUMERICAL RESULTS AND DISCUSSION 

We have solved numerically eqs. (p^)^jl^), the results 
may be summarized as follows. In Fig. [ijphase diagram 
as a function of n and V for different values of U for the 
two-dimensional square lattice at zero temperature. The 



FIG. 2. (T-V)-phase diagram at quarter filling for several 
values of U : U = 1, 2 and 3 corresponding the dotted, solid 
and dashed lines, respectively. The inset shows (U-V)-phase 
diagram for n — 1/2 and T = 0. 



half-bandwith W was taken as unit of energy (for the 
square lattice z = 4 and W = it). Due to the electron- 
hole symmetry we consider only < n < 1. From Fig. |l| 
one can see that in the two regions (n <£L 1, n < 1) the 
influence of U on the boundary between the charge or- 
dered and the homogenous disordered phase is different: 
away from half filling (n < n* ss 0.67) the on-site inter- 
action U has little effect on the critical value V c , while 
for Til < n < 1 the critical V c strongly depends on U. In 
Ref.Ej Dho et al found that CO exists over a broad dop- 
ing range (0.44 < x < 0.8) in La 2 -2a;Sri + 2 2: Mn207. It 
is worthwhile to note that for large U in the same filling 
region the CPA values of V c are close the minimum value 
(V c has the minimum at n* w 0.67). Although the mech- 
anism of the charge ordering in the layered manganitcs 
is more complex than one induced by a nearest neighbor 
Coulomb repulsion, from the above result we may specu- 
late that the CPA is able to describe the CO boundary in 
this compound. In addition, the advantage of the CPA is 
that by using this simple approach one can easily obtain 
phase diagrams in the (V^-n)-plane for arbitrary U as well 
as in the ([/-y)-plane at arbitrary n. As an illustration 
of our approach in the following we consider CO tran- 
sition on the square lattice at quarter filling (n = 1/2). 
The inset in Fig. || shows the (U-V)-phase diagram at 
zero temperature. We compare CPA result to the one 
obtained by other methods. At U — oo the critical .walue 
V c = 0.195IF and V c = 0.172W obtained in Ref.E3 by 
the slave boson approach with a constant and the actual 
DOS po(w) respectively, is in a good agreement with our 
result V c = 0J218W. At U = 2W V c = 0.66W was ob- 
tained in RefH for the EHM in infinite dimension by the 
numerical renormalization group (NRG) method, while 
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FIG. 3. (T-V)-phase diagram for U = 2 at various values 
of n : n = 0.3, 0.5, 0.65 and 0.8 corresponding to the solid, 
long-dashed, dotted and dashed lines, respectively. The inset 
shows the lattice occupancies ua (solid line) and ns (dashed 
line) as a function of V for U — oo and T = 0. 



FIG. 4. CPA results at quarter filling for the A and B spec- 
tral functions p(u>) (solid and dashed lines, respectively) for 
U = co, T — and several values of V. Below the criti- 
cal V c ~ 0.218, the spectral functions of both sublattices are 
equal. 



the CPA-result for our 2D lattice is V c = 0.237W. 

In Ref.lI3 McKenzie et al argue that in the large U limit 
at quarter filling the charge ordered phase is destroyed 
below a critical non-zero value V c , of the order of t, while 
we can now show in the inset of Fig. || that the critical 
value V c is almost independent of U, therefore V c is of 
the order of t for all values of U. Fig. |] shows the (T-V)- 
phase diagram for different values of U. For each value of 
U reentrant behavior as function of temperature is seen 
for some region of V. For V in this range, the ground 
state is homogeneous but a charge ordered phase exists 
at intermediate temperature. Note that reentrant transi- 
tion apparently is not obtained by Hartree-Fock approxi- 
mation in the EHM without electron-phonon interaction. 
Fig. U shows the (T-y)-phase diagram for U — 2 at dif- 
ferent band fillings n. The reentrant behavior is seen 
clearly for the values of n where V C (T = 0) is larger (e.g. 
n = 0.3, n = 0.8), while for values of n where V C (T = 0) 
is small (e.g. n = 0.65), the reentrant behavior is not 
seen. 

In the inset in Fig. || we present the CPA result for V 
dependence of the sublattice occupancies ua and ns and 
hence the charge order parameter (ua — Ub) in the strong 
correlated limit U — > oo at zero temperature. The tran- 
sition is clearly continuous in contrast to the result for 
U = 2, T = in Ref.ll where the NRG method gives 
a first order phase transition. Note that in the strong 
correlated limit U — > oo from eq. (11) one can find an 
analytic expression for G(x,oj). Then eqs. ( |l2| ) have to 
solved self-consistently to find ha, ns and /i. For a given 
set of these parameters the total energy of the system for 
various states can be calculated. The state with the low- 



est energy is the true ground state and determines the 
spectrum. The CPA result for A and B spectral func- 
tions for U — oo, T — are shown in Fig. ||. When 
V < V c ~ 0.218 the A and B spectral functions are 
identical and they are independent of V by shifting all 
the one-electron energy levels and the chemical poten- 
tial by 2V. For V > V c the A and B spectral func- 
tions change and each spectrum splits in the upper and 
lower subbands. On increasing V the weight of the lower 
subband in the A-spectrum increases, while the one in 
the B-spectrum decreases; the subbands in the spectrum 
become narrower due to the reduced hopping of elec- 
trons in the charge ordered phase. A charge ordering 
gap opens and it is given by A — 2y^l6V 2 x 2 + 7 , where 

7 = " f ) " ifX 1 " ,J f)l In the limit of 

large V, perfect charge order evolves: ua — ► 1, — » 0. 
Therefore in this limit x — > 1/2 and the gap between two 
peaks A is given in CPA by JlV, as compared to 2V ob- 
tained by the DMFT in RefJi 



IV. CONCLUSIONS 

In this paper we have applied the CPA to study charge 
ordering in the extended Hubbard model. Within this 
approximation one can obtain the critical value V c as a 
function of temperatute T, on-site Coulomb repulsion U 
and band filling n. To examine the CPA results we con- 
sider the charge ordering transition on the 2D square lat- 
tice at quarter filling using a semi-elliptical DOS. It was 
shown that for all values of U the charge ordered phase is 
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destroyed below a critical non-zem, value V c of the order 
of t. Like previous results in Ref.ull3, we find a parameter 
region where the model shows reentrant behavior. The 
reentrant transition is also observed at other band filling, 
as it was. experimentally found in the layered manganites 
in Ref.uJ. In the strong correlation limit U — * oo at zero 
temperature the CPA gives a continuous transition. 
Now the CPA is known to give good results of one- 
particle properties in a wide range of systems. To study 
the CO boundary phase in the EHM the CPA has the 
advantage over DMFT of being analytically simple, and 
over Hartree-Fock approximation (small U) and slave bo- 
son approach (U — ► oo) of being able to describe the 
whole range for the on-site interaction U. Of course in 
common with the alloy CPA the imaginary part of the 
self-energy does not vanish at the Fermi level at T = 0, 
so we do not obtain a true Fermi liquid and we only ex- 
pect our CPA to study CO in the EHM well at finite 
temperature. 

The calculation presented here can also be applied to 
the lattice of higher dimensions, or to the EHM in the 
presence of a weak magnetic field. To include magnetic 
phases and cluster effects one has to go beyond the usual 
alloy CPA. This is left for our future work. 
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